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Based on a Fermi liquid picture, the temperature effect on the impurity-induced spatial modulation 
of local density of states (LDOS) is investigated for the d-wave superconductor Bi2Sr2CaCu208, in 
the context of scanning tunneling microscopy (STM) . It is found that stripe-like structure exists even 
in the normal state due to a local-nesting mechanism, which is different from the octet scattering 
mechanism proposed by McElroy et al. [Nature 422, 592 (2003)] in the d-wave superconducting 
(dSC) state. The normal-state spectra, when Fourier-transformed into the reciprocal space, can 
reveal the information of the entire Fermi surface at a single measuring bias, in contrast to the 
point- wise tracing proposed by McElroy et al. This may serve as another way to check the reality 
of Landau quasiparticles in the normal state. We have also re-visited the spectra in the dSC state 
and pointed out that, due to the Umklapp symmetry of the lattice, there should exist additional 
peaks in the reciprocal space, but experimentally yet to be found. 
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I. INTRODUCTION 

A fundamental question on the high- Tc cuprates re- 
tnains after more than fifteen years of the discovery of 
the material. It is still not clear whether the cuprates 
are systems of Fermi liquid (FL), non-FL with exotic or- 
ders such as the stripes [1], or systems with more in- 
tricate co-existence of different states of matter [2]. In 
fact, this question arises in both their superconducting 
and normal states. Besides, it is not unusual that differ- 
d ent experimental probes give different implications - the 
experimental findings are not yet converged. 

A substantial progress in the STM measurement has 
made it another route to this problem. Not only 
STM looks directly into the real space, but also it can 
be readily connected to the reciprocal space, so-called 
the Fourier transformed STM (FT-STM). Data of the 
low (and fixed) temperature STM on Bi2Sr2CaCu208 
K> (BSCCO) was claimed to be an excellent manifestation 
^ of the FL behavior [3,4]. But, since the cuprates are such 
5^ involved systems, one should be more careful to nail down 
the conclusion. Whether the observed STM modulation 
[5] is solely the Friedel stripe arising from the quasipar- 
ticle (QP) interference, or the Zaanen-Kivelson stripe [1] 
coexisting with the Friedel stripe, is actually an issue 
still in debate (see the contradictory data of Refs. [6] and 
[7]) [8]. Even if the stripes can be attributed to the QP 
interference alone, it is still crucial to ask how well do 
the quasiparticles behave? How is the extension of the 
"Fermi-arc" , and up to what temperatures they survive 
[9]? So far there has been few finite-temperature and 
normal-state STM studies on the cuprates [10]. These 
are the major concerns of the present paper. 

Based on the FL scenario, an abstract model named 



the "octet" scattering model [5] has successfully as- 
cribed the experimentally observed FT-STM peaks of the 
LDOS modulation to the quantum interference of the 
QPs. Later concrete single-impurity scattering calcula- 
tions [11,3,4] also supported that. The occurrence of FT- 
STM peaks and their evolution with the bias change give 
information of the Fermi surface (FS) of the measured 
system which are consistent with previous results from 
the angle-resolved photoemission spectroscopy. This pro- 
vides a strong identification of a Friedel stripe out of 
the Zaanen-Kivelson stripe. But nevertheless, some weak 
non-dispersive peaks, which could be due to a coexisting 
Zaanen-Kivelson stripe, may also exist [6~8]. 

It is important to investigate how far the FL picture 
can be pushed, especially when the temperature is raised 
and the system enters the normal state. We thus have 
performed similar FL-based calculations at different tem- 
peratures to study the Friedel stripes, for the reference 
of future experiments in checking the validity of the FL 
picture. It is found that stripe-like structure exists even 
in the normal state. This may be counterintuitive to 
the octet model, since the octets should vanish in the 
normal state. We argue that apart from the octets, a 
local-nesting property of the FS (especially important in 
the normal state) can also give rise to sizable joint-DOS, 
and hence sizable QP scattering. 

In Sec. II, we formulate the problem as the scattering 
from a single impurity. In Sec. Ill, we present the calcu- 
lated spectra at different temperatures and the essential 
features are remarked. In Sec. IV, the physical origin 
of the normal-state spectra are illucidated. It is noted 
that in the normal state, the information of the entire 
FS can be revealed from the data of a single measuring 
bias. This is in contrast to the point- wise tracing out in 
the dSG state. Sec. V is a closing section. In Appendix, 



we make a comment to the present FT-STM data, point- 
ing out that some additional peaks should exist in the 
measurement. 



II. FORMALISM 

Since the single-impurity scattering model was proved 
to be an excellent start for understanding the STM fea- 
tures in the dSC state, we proceed to study the finite- 
temperature phenomena based on it. We will be inter- 
ested in those regions away from the impurity neighbor- 
hood. We consider the following Hamiltonian 

H = Ho + Hi, (1) 

where Hq is the usual BCS mean-field Hamiltonian 



Ho = ^^k4^cka + [^^(i^^-ki + H-c- 



(2) 



k,(7 



with and Ak the band dispersion relative to chemical 

potential and the superconducting gap function respec- 
tively, and Hi is the part associated with an impurity at 
site 0, 



H 



(3) 



Parameter 6t is the deviation of local hopping, SA is the 
deviation of local pairing potential, and Vq is the on- 
site impurity potential. This is a simplified model but 
nevertheless should be enough to bring out the essential 
features of interest. 

It is convenient to apply the Nambu representation for 
the reciprocal and real-space operators: 



,t 



and 



(4) 



which are related to each other via the Fourier transfor- 
mation 



Ck 



(5) 



Here the sum is over all the lattice sites r j . Define a 2 x 2 



energy matrix 



thus 



Similarly 
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At 



Ak 
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Ho = YclikCk. 



(6) 



(7) 



Hi = Yc}uijCj+R.c, 



(8) 



where the matrix elements Uij are to be given later. 

Since we are interested in the real-space STM spectra 
which measure the LDOS, we need to know the real- 
space, cqual-sitc, singlc-particlc Green's functions. De- 
fine a single-particle Green's function matrix 



G{ri,Tj,T) = -{TrCi{T)C}{0)), 

the LDOS is then given by 
1 

- G22(r, r, -iu, 



(9) 



-Im[G'ii(r, r, iuj„ 



-^0+)], 



(10) 



where Ga/3 is an element of the 2x2 Green's function ma- 
trix in (9), being Fourier-transformed to the Matsubara- 
frequency space 



Jo 



(11) 



Following the standard technique, the full G in (9) can 
be expanded in terms of Hi given by (8), 

G{v,.Vj,T)=G\TuTj,T) 

+ ! dTl{TrCi{T)Y,Cl{Tl)uuCi{Tl)C]m 



+0{u^) 



(12) 



Vj,T) is the mean- field 



where. G"(rj,rj,T) = G'^(r: 
"non-interacting" G when Hi = Q. 

The terms included in Hi in (3) are given explicitly 
here. For BSCCO, wc consider a square lattice of lattice 
constant a. In the case of a single, extended and weak 
impurity, we consider local deviations 5t\ and (5Ai which 
couple the impurity site and its nearest neighbors, and 
5t2 and 5/S.2 which couple the impurity's nearest neigh- 
bors and its next nearest neighbors. Consequently, there 
are 17 non- vanishing u matrices in (8): 



^0,0 







-Vo 

5ti -SAi 
SAi -Sti 



Uo,±ax = 



5ti 5A-1 
-SAi -Sti 



^ — ax, — 2ax — ^ay,2ay — ay,— 2ay 

5t2 SA2 ' 
-SA2 -St2 

'^ax,ax:kay ~ '^—ax,—ax^ay ~ "^ay^ayzkax ~ '^—ay,—ay:kax 



5t2 

6A2 



-SA2 

-6t2 



(13) 



Hamiltonians similar to this have been successfully used 
by Tang and Flatte [11] to explain the resonant STM 



spectra for Ni doped BSCCO, and by Wang and Lee [3] 
and Zhang and Ting [4] to explain the energy-dependent 
modulation of the FT-STM spectra on superconducting 
BSCCO. Since we consider the weak impurity scattering 
limit, G in (12) is only calculated up to the first order 
of u (the Born limit). The first order term has already 
included the essential interference effect of the; QPs. A 
strong impurity is expected to give new features (such as 
a bound state) only at the immediate neighborhood of 
the impurity. Eq. (12) is then reduced to 

G{vi, rj,iun) = G°(ri - r^, iw„) 

+ ^G°(rj - rk,iu)n)uk(G°{re - rj,ia;„). (14) 

The first term on the right is translationally invariant. 
Spatial variation of the LDOS at a fixed energy comes 
only from the second term. It is related to the STM 
spectrum discussed in Rcf. [5] via = Vj = r and icun 
eV + iO~^, where e = |e| is the electron charge and V is 
the bias voltage. Its Fourier transform is the FT-STM 
spectrum. We will be discussing the second term on the 
right throughout this paper. 

Note that the thermal Fermi distribution function 
never appears in the full Green's function G here. It 
arises only when the Matsubara-frequency sum is in- 
volved, i.e., when the impurity is dynamic and inelastic 
scattering occurs. Since we have assumed an elastic im- 
purity, our Green's function depends on temperature only 
through the gap magnitude A(T) which should change 
with temperature. 

In our calculation, we have used a 800 x 800 square lat- 
tice with the impurity at the center. We have chosen a 
simple but reasonable impurity potential, 26ti = 45t2 = 
—26Ai = — 4i5A2 = Vq and have assumed that these 
scales are small and in the perturbative limit. For ^k, 
we use a tight-binding band, = ti{coskx + cosky)/2 + 
t2CoskxCosky + t3{cos2kx + cos2ky)/2 + t4{cos2kxCosky + 
coskxCos2ky)/2 + t5Cos2kxCos2ky — fi (lattice constant 
a = 1), with ti_5 = -0.60,0.16,-0.05,-0.11,0.05 eV 
and chemical potential ji ~ —0.12 eV, appropriate for an 
optimally-doped BSCCO [12]. In addition, the supercon- 
ducting gap is taken to be Ak = A(T)(cos kx — cos ky)/2 
with A(T) the temperature-dependent gap magnitude. 
Besides, we have introduced a finite broadening 7 = 2 
meV to the Green's function, such that eV + i0+ is re- 
placed by eV + ij. 

III. STM AND FT-STM SPECTRA 

In Fig. 1, we present the temperature evolution of the 
real-space STM spectra at two different negative bias 
voltages. The case of positive bias voltage will not be 
discussed as they are qualitatively the same. The gap 
magnitudes are taken from A(T = 0) = 44 to A(rc) = 
meV to simulate the transition from the superconducting 



to the normal state. At a distance of several lattice con- 
stants away from the impurity, oscillating Friedel stripes 
are seen in all panels, even in the case of normal states 
[A(r) = 0]. Spectra with a similar ratio of e\V\/A{T) 
share a similar behavior, such as those in Fig. 1(b), (c) 
with e\V\/A{T) = 25/44,15/26 - 0.57. On the other 
hand, spectra at zero A(T) is robust at the change of 
V [see Fig. 1(g) and (h)]. Comparing the relative in- 
tensities of the modulations, we see that the strongest 
modulations (at a fixed bias voltage) appear at temper- 
atures near Tc. One more important feature to note is 
that the different ripples live in well-separated patches of 
space. The ripples do not have large- area overlaps. 

The Fourier-transforms of the above spectra are given 
in Fig. 2. Roughly speaking, there are two regimes for 
the spectra in the A(T) ^ superconducting state, as 
classified by the ratio e\V\/A{T). When e|y|/A(r) < 1, 
there are local peaky structures that have their locations 
describable by the octet model [see Figs. 2(a)-(c)], in 
agreement with previous studies [3,4]. Previous stud- 
ies discussed the cases with different eV and a fixed 
A(T) A(0), while in our case A(T) is varied (by chang- 
ing the temperature). In comparison of Fig. 2(b) and (c) 
that have a similar ratio e|F|/A(T) ^ 0.57, one again 
sees that the essential features of the spectra depend 
mainly on the ratio e|F|/A(T), as noted before in the dis- 
cussion of the real- space spectra. In Fig. 3, a closer look 
of the locations of the interference peaks in Figs. 2(a)- 
(c) is given. When e\V\/A{T) > 1 [see Figs. 2(d)-(f)], 
some extra peaky structures beyond the description of 
the octet model appear near q = 0. These new struc- 
tures are not expected to be observable in real BSCCO 
compound because they are related to the maximum gap 
part of the STM spectra (see discussion in Sec. IV), which 
in turn are highly inhomogeneous in space. As A(T) is 
further decreased, the strong peak at q = is suppressed 
and vanishes at entering the normal state. 

In the normal state, the FT-STM spectra are reduced 
to some neat ridges instead of peaks [see Figs. 2(g) and 
(h)], and they are rather robust against the change of 
bias. They should be readily observable in practice. 

As the interference peaks in the dSC states (in the 
regime of e\V\/A{T) < 1) were well documented in the 
literature [3,4], we will only give a supplementary com- 
ment on it in Appendix. In Sec. IV we pay special atten- 
tion to the normal-state spectra. 

IV. ORIGIN OF THE NORMAL-STATE 
SPECTRA 

In this section, we show that the normal-state FT-STM 
spectra have an intimate relationship with the underlying 
FS. 

The pronounced feature of the FT-STM in the dSC 
state is that the peaks have locations depending on the 
bias. As long as the ratio e|y|/A(T) is small enough 



to stay away from the maximum gap region, locations 
of the peaks are more or less as described by the octet 
model [3] - a model that assumes the dominant QP scat- 
tering comes between regions of high DOS on the FS 
(the octets). These are regions with the smallest velocity, 
which appear at the tips of the banana-shape constant 
energy contour (|^kP + 1 Ak(T)|^ = constant). The agree- 
ment of the low-temperature experimental observation 
with the picture was claimed to be good [5]. Later spe- 
cific single-impurity scattering calculations [3,4] (similar 
to the one done in Sec. II) also supported this picture. In 
the normal state, such large DOS octets no longer exist. 
However, we point out that there exists a different mech- 
anism which can also cause a substantial joint-DOS, and 
leads to distinguished structures in the q space. 

The upper panel of Fig. 4 shows a typical normal-state 
FT-STM spectrum in an extended Brillouin zone (BZ). 
The underlying FS is also shown in the lower panel. Com- 
paring the two panels, it is readily seen that the "ridges" 
in the spectrum are of the same shape as the FS, but 
having twice the size, and differently oriented branches 
overlap together. The occurrence of the ridges can be un- 
derstood from the scattering wave vectors drawn in the 
lower panel. Those wave vectors, which are pivoted at 
q = {rmr/a,mr/a), n,m & integer, are special in the sense 
that they joint locally parallel segments of the FS, i.e., 
they possess a weak "local-nesting" [13] (such nesting is 
the weakest possible type of nesting, or it is "marginal"). 
As a result, stripes in the real-space are still understood 
as due to the scattering on the FS. 

In the normal state, there is no small and well- 
separated large DOS regions, therefore the stripes in real- 
space are smoothly deformed unlike those in the dSC 
state. An implication from the above understanding to 
the dSC state is, the new structures at higher eV/A{T) 
in Fig. 2 (which were not seen by the experiment) come 
actually from the maximum gap regions. In either state, 
quasiparticles come from different directions Vi;- to hit the 
impurity, get bounced back, and get interfered with some 
transition wave vector q. Since and q are not simply 
related, the orientations of the interference ripples in the 
patches of space are also not simply related to their orien- 
tations from the impurity. It is also obvious that within 
the same patch of space there will be no two crossing rip- 
ples. Therefore, to explain the "checkerboard" pattern 
(overlapping ripples) seen in the experiments using this 
picture, one needs to consider the existence of a dilute 
concentration of impurities in the system [14,15]. 

V. CONCLUDING REMARKS 

Our perturbative approach should not be a concern re- 
garding the validity of our discussion. In this paper, it 
is not our purpose to investigate the local bound states 
at the vicinity of the impurity. We are interested only 
in those regions remote from the impurity, where the in- 



terference between scattered quasiparticles is expected to 
be the dominant process. The approach has included the 
interference effect lucidly. 

In the literature, it was proposed that temperature 
could be a detrimental factor to the quasiparticles [9]. 
If this is true, the normal-state stripe described in this 
paper is expected to vanish. Stripes of a different origin 
may exist [1] , but they are likely to differ in many aspects, 
such as the orientation or temperature dependence of the 
stripes. In our case, the real-space normal-state stripes 
[14] and the corresponding ridges in the FT-STM spec- 
tra reflect the information of the entire Fermi surface. 
If there are destructions of quasiparticle states at some 
Fermi surface segments, the corresponding segments of 
the ridges in the FT-STM spectrum should also be de- 
stroyed. The presence of a pseudogap in the normal state 
should also be signified as some missing segments of the 
ridges. Our study interpolates the dSC and normal-state 
STM, using an approach supported by the experiment in 
the dSC state. The "fingerprints" of the Fermi liquid un- 
der the STM probe, in the course of the superconducting- 
to-normal transition, are enumerated for future experi- 
ments. 
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APPENDIX: UMKLAPP SYMMETRY 

In this Appendix, a brief comment is given to the ex- 
perimental data in Ref. [5]. If one considers the octet 
model in the dSC state more carefully, it will be interest- 
ing to note the existence of additional peaks, which are 
not yet realized in current experimental data. 

In Fig. 5, for a certain bias voltage, we have illus- 
tratively shown a few direct and Umklapp scattering 
wavevectors in the octet model [16]. In addition to the 
scattering vectors within the first Brillouin zone (direct), 
there are also vectors connecting octets in different Bril- 
louin zones (Umklapp). When all these vectors are taken 
into account, there should exist FT-STM peaks as shown 
in the lower panel of the figure. Take an example, a new 
peak at the vector q'^ is seen along the qy axis (when 
the bias voltage exceeds some value). Peak at gs was 
reported by Ref. [5], but the corresponding peak was 
not. The more disturbing matter is, there also exists a q'^ 
peak which may come close to ^2,6 peaks (or vice versa 
speaking). 

It may seem at odds that a mathematical Fourier trans- 
form of experimental data on a discrete lattice does not 



automatically possess the Umklapp symmetry. The rea- 
son may be duo to the fact that practical STM actually 
probes a space more continuous than the underlying lat- 
tice. It scans in steps smaller than the lattice constants, 
and identifies the locations of individual atoms only later, 
by looking at the modulation profile. The probe tip does 
not jump from one atom to another right on-top! 



FIG. 1. Temperature dependence of the LDOS modulation 
surrounding the impurity. All panels axe plotted within a 
square region, rx,ry = —20 ~ 20a, with the impurity at the 
center. The left and right columns show the spectra at bias V 
= -15 and -25 mV respectively, and from top to bottom the 
gap magnitude A(T) = 44, 26, 10, and meV respectively. 
For better visualization, spectra in different panels are plotted 
in different intensity windows from —x (deepest blue) to +x 
(deepest red) as specified by the color scale at the bottom. 
The individual value of x is given at the lower-right corner of 
the individual panel. Those regions of the deepest blue/red 
have intensities below/above the window. 
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FIG. 2. Temperature dependence of the FT-STM spectra 
(corresponding to the real-space results in Fig. 1) in the first 
Brillouin zone qx,qy ~ —it /a ~ Tv/a. Left and right columns 
are for V = —15 and -25 mV respectively. From top to bot- 
tom, A(T) = 44, 26, 10, and meV respectively. The spec- 
trum in each panel is plotted in an intensity window from xi 
to Xu- The values of xi and specified at the bottom of 

each panel, and the color scales are shown at the bottom of 
the figure. 



FIG. 3. This figure shows that the positions of the inter- 
ference peaks (indicated by arrows) depend mainly on the 
ratio e|F|/A(T). Spectra in Figs. 2(a)-(c) are scanned along 
q/a = (0.567r, 0.567r) -> (0,0) -> (0.87r,0). It is seen that 
Figs. 2(b) and (c), which have a similar e|V^|/A(T) ratio, have 
similar peak positions and intensities. For easier comparison, 
the intensity of the spectrum in Fig. 2(a) is enlarged by a 
factor of 2. 



FIG. 4. (Upper panel) The FT-STM spectrum of a nor- 
mal-state system at a bias of -23 mV. The ridge-like structure 

is a typical feature of all normal-state spectra. (Lower panel) 
Several of the wavevoctors on the above ridges are shown in 
the extended Brillouin zones. They are wavevectors having a 
local- nesting property, i.e., they joint locally parallel parts of 
the Fermi surface (as illustrated by pairs of red lines). 



FIG. 5. (a) For illustration, a few of the scattering wave 
vectors in the octet model (the octets are indicated as cir- 
cles) are shown. Wavevectors discussed by McElroy et al. 
[5] are shown in blue, while the Umklapp wavevectors (which 
were not discussed in Ref. [5]) are in red. (b) Locations of 
all the wavevectors connecting all the octets. The Umklapp 
wavevectors are in red. 
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